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CN . Abstract 

["t"! ■ This paper proposes to extend the discrete Verhulst power equilibrium approach, previously suggested in 

Qj, to the power-rate optimal allocation problem. Multirate users associated to different types of traffic are 
Q , aggregated to distinct user' classes, with the assurance of minimum rate allocation per user and QoS. Herein, 

Verhulst power allocation algorithm was adapted to the single-input-single-output DS/CDMA jointly power-rate 
^ ■ control problem. The analysis was carried out taking into account the convergence time, quality of solution, in 

' terms of the normalized squared error (NSE), when compared with the analytical solution based on interference 
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, allocation methodology. 
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matrix inverse, and computational complexity. Numerical results demonstrate the validity of the proposed resource 



Index Terms 

Power-rate allocation control; SISO multirate DS/CDMA; discrete Verhulst equilibrium equation; QoS. 

I. Introduction 

In the last years many efforts has been spent trying to find the best resource allocation algorithm that 
could be easy applied to DS/CDMA communications systems. The Foschini and Miljanic [|2|| studies 
can be considered as a foundation of many well-known distributed power control algorithms (DPCA) in 
scientific literature, because they try to solve an ordinary differential equation (ODE, eq. (1) in |l2l), which 
with some minor alterations is also considered in many other subsequent studies. Therefore, a new ODE 
can lead to a new algorithm, more promising in several aspects, such as, convergence, proximity to the 
optimum value, and sensibility to estimation errors as well. 

With this context in mind, the work in |[T1 proposed and analysed a new ODE for the DPCA based 
on the Verhulst equation |l3l. The discrete version of the Verhulst population model is more diffused in 
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the literature and it is called logistic map. The logistic map was studied thoroughly by R.M. May in dUl. 
Von Neumann and Ulam (Si also studied the logistic map and they evaluated the possibility to use it as 
random generator number, which is gotten in certain conditions. 

The Verhulst model was initially designed to describe population growth of biological species with food 
and physical space restriction. With the successfully mathematical model adaptation to power control in a 
single-rate DS/CDMA systems proposed in [|T]|, this work suggests an expansion of the Verhulst approach 
to other optimization problems, such as the power-rate allocation problem, and its multi-objective versions, 
due to DPCA performance x complexity gain when compared to the classical algorithms such as Foschini 
or sigmoidal. 

In this paper we have adapted the Verhulst approach to the power-rate allocation problem with multirate 
QoS associated to different types of traffic (based on user classes), and satisfying the minimum rate 
allocation per user requirement. Hence, the Verhulst power allocation algorithm of |[T1 was adapted to the 
power-multirate control problem. The analysis was carried out taking into account the convergence time, 
quality of solution when compared with the analytical solution based on interference matrix inverse, and 
computational complexity. 

This paper is organized as follows: Section |Il] gives an overview of the power control classical solution 
and how it is adapted to multirate problem. In Section Hn] the power-rate Verhulst algorithm is proposed. 
Numerical results with corresponding simulation parameters setup are treated in Section |IVl Finally, the 
conclusions are offered in Section |Vl 

II. Power and Rate Allocation Problem 

In a multiple access system, such as direct sequence code division multiple access (DS/CDMA), the 
power control problem is of great importance in order to achieve relevant system capacity and throughput. 
The power control problem can be solved by a vector that contain the minimum power to be assigned in 
the next time slot to each active user, in order to achieve the minimum quality of service (QoS) through 
the minimum carrier to interference ratio (CIR). 

In multirate multiple access wireless communications systems the bit error rate (BER) is often used as 
a QoS measure and, since the BER is directly linked to the signal to interference plus noise ratio (SNIR), 
we are able to use the SNIR parameter as QoS measurement. Hence, associating the SNIR to the CIR at 
time slot n results: 

6,[n] = -^xr,[n], n = 0,l,...N (1) 

Ri[n\ 

where 6i[n] is the SNIR of user i at the nth iteration, Rc is the chip rate, Ri[n] is the data rate for user 
i, Ti[n] is the CIR for user i at iteration n, and is the maximal number of iterations. From ([T]) we are 
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able to calculate the data rate for user i at iteration n: 

Rr 



Rih] 



SAn] 



n 



The CIR for the ith user can be calculated as [m, 

Pi[n]gii[n] 



r^n] 



K ' 

E Pi[n]gij[n] + cr2 



0,1,. ..,iV 



(2) 



(3) 



where Pi[n] is the power allocated to the ith user at time slot n and is bounded by [Pmin; ^max], the 
channel gain (including path loss, fading and shadowing effects) between user j and user (or base station) 
i is identified by gij, K is the number of active users in the system, and af = = is the average 
power of the additive white Gaussian noise (AWGN) at the input of ith receiver, admitted identical for 
all users. Therefore, in DS/CDMA multirate systems the CIR relation to achieve the minimum rate can 
be calculated to each user class as follows BTJ: 

I'm in ^ 



Rr 



i = l---L 



(4) 



where F^j^^ and i^^j^ is the minimum CIR and minimum user rate associated to the £th user class, 
respectively, 6* is the minimum (or target) signal to noise ratio (SNR) to achieve minimum acceptable 
BER (or QoS), and L is the total number of user classes in the system (voice, data, video, and so on). 
Besides, the power allocated to the kth user belonging to the ith class at nth iteration is: 



k=l---Kf- 



(5) 



hence, the total number of active users in the system is given by K = Ki U . . . U U . . . U Ki^. Note that 
indexes associated to the K users are obtained by concatenation of ascending rates from different user's 
classes. Hence, Ki identifies the lowest user's rate class, and the highest. 

The K X K channel gain matrix, considering path loss, shadowing and fading effects, between user j 
and user i (or base station) is given by: 



G 



gu gi2 
g2i g22 



giK 
g2K 



gxi gK2 ■ ■ ■ gxK 

which could be assumed static or even dynamically changing over the optimization window (A^ time 
slots). 
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Assuming multirate user classes we are able to adapt the classical power control problem to achieve 
the minimum rates for each user, simply using the Shannon capacity relation between minimum CIR and 
minimum rate in each user class, resulting: 

r^, = 2<." - 1 (6) 

Now, considering a x K interference matrix B 



0, i = j 



Bj,- = < r- ■ a - (7) 

where Fj can be obtained from dH), taking into account each rate class requirement, and the following 
column vector x 1: 

u. = (8) 

9ii 

we can obtain the analytical optimal power vector allocation simply by matrix inversion as: 

p* = (I-B)U (9) 

if and only if the maximum eigenvalue of B is smaller than 1 (81; I is the K x K identity matrix. In this 
situation, the power control problem shows a feasible solution. 

Herein, the classical power allocation problem is extended to incorporate multirate criterium in order 
to guarantee the minimum data rate per user class. Mathematically, we want to solve the following 
optimization problem: 

min p = [pI . ..p]^^, . . . . . -Pk,, ■■ -Pki] 
s.t. Pin < pi < PL. (10) 
= R^.^^ yk e Ki, and V£ = 1, 2, ■ ■ ■ L 

III. Verhulst Power-Rate Optimization Approach 

The Verhulst mathematical model was first idealized to describe population dynamics based on food 
and space limitation. In fll that model was adapted to single-rate DS/CDMA distributed power control 
using a discrete iterative convergent equation as follows: 

. St 

where pj[ri + 1] is the user i power at the n + 1 iteration, a E (0; 1] is the Verhulst convergence factor, Si[n\ 
is the ith user' SNIR at iteration n, 5t is the minimum SNR for the iih. user that guarantee a minimum 
QoS in terms of performance (BER). 



Pi[n +!] = (! + a)pi[n] - a 



1,---,K (11) 
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The recursion (fTTl) can be effectively implemented in the ith mobile unit since all necessary parameters 
a, the QoS level given by delta*, the transmitted power Pi[n], except 6i[n\, can be considered known in 
the mobile unit i. The SINR 6i[n] can be obtained only at the correspondent base station that demodulates 
the signal from user i. In this way, the BS estimates Si[n], quantizes it in a convenient number of bits, 
and transmits this information to the ith user through the direct channel. Thus, (fTTl) depends on local 
parameters just allowing that the power control works in a distributed manner, i.e., each one of the K 
links (mobile terminals to base station) carries out separately the respective power control mechanism, 
justifying the name distributed power control algorithm (DPCA). 

Equation (fTTl) gives a recursive power update, close to the optimal power solution after iteration. 
However, originally it does not consider the rate requirements in a multirate environment. In order to 
achieve the QoS to each user class, (fTT|) must be adapted to reach the equilibrium lim Pi[n] = p* when 

n— >oo 

the power allocated to each user satisfies the minimum rate constraint given in (fTOl) . Hence, the recursive 
equation must be rewritten considering SNIR per user class, via (HI), in order to incorporate multirate 
scenario. The minimum CIR per user class is obtained directly by dH). In this way, considering the 
relation between CIR and SNIR in a multirate DS/CDMA context, we propose the equation below in 
order to iteratively solve optimization problem in (fTOl) : 



6,[n] = FxT,[n] (12) 

= Fx— ILJM^ ,i=l,...,K 

E PM9ij[n] + cr2 



where F is the spreading factor per user class, given by: 



(13) 



-^min 

Note that the CIR of ith user at the nth iteration is weighted by spreading factor; so the corresponding 
SNIR is inversely proportional to the actual (minimum) rate of the ith user of the ^th class. 

A. Quality of Solution x Convergence Speed 

The quality of solution achieved by iterative Verhulst equation (fTT|) is measured by how close to the 
optimum solution is p[n], and can be quantified by means of the normalized squared error (NSE) when 
equilibrium is reached. The NSE definition is given by: 



NSE[n] = E 



IpN -P*f 

II *l|2 

Hp II 



where || ■ |p denotes the squared Euclidean distance to the origin, and E[-] the expectation operator. 
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On the other hand, the convergence speed in Verhulst equation is dictated by the parameter a. Hence, 
for small values of convergence factor, i.e., a — i- 0, the convergence is slow, but the NSE is very small 
after iteration, when compared with the opposite configuration: the convergence is fast when a — )• 1, 
but the NSE is a concern. So, in order to accelerate convergence, we propose two adaptive criteria for 
convergence factor a when iterations evolve, based on a) SNIR to target SNR difference, and b) tanh 
mapping for this difference, as following: 

a) ai[n] = min ja^ax; ^'^'^^ ^} — + aminj , (14) 

b) ai[n] = maxlttinin; tanh(|(5i[n - 1] - 6*\)} , (15) 

with Omin = 0.1 and Omax = 0.95. 

IV. Numerical Results 

Simulations were carried out through the MatLab ver.7.3 platform, with system parameters indicated in 
Table m For all simulation results discussed in this section, it was assumed a retangular multicell geometry 
with a number of base station equal to 4 and mobile terminals uniformly distributed. A typical placement 
for mobile terminals (mt) and base stations(BS) is provided in Fig. \T\ Besides, the rate assignment for 
all multirate users was considered uniformly distributed as three submultiple rates of chip rate, -Rmin = 
[128' 32' Tel-^c [bps]- 

A number of mobile terminals ranging from = 5 to 30 was considered, which experiment slow 
fading channels, i.e., the following relation is always satisfied: 

Tsiot < (At)^ (16) 

where Tsiot is the time slot duration, and (At)^ is the coherence time of the channel^. This condition is 
part of the SINR estimation process, and it implies that each power updating accomplished by the DPCA 
happens with rate of T^^t^ assumed here equal to 1500 updated per second. The recursion in (fTTj) should 
converge to the optimum point before each channel gain gij experiments significant changing. Note that 
satisfying (fT6l) the gain matrices remain static during one convergence process interval. 

In all simulations the entries values for the QoS targets were fixed in 5* = 4 dB, the adopted receiver 
noise power for all users was F„ = —63 dBm, and the gain matrix G entries had intermediate values 
between those used in flU and [|71. Furthermore, in order to evaluate and highlight the different aspects 
and features of the proposed resource allocation methodology, the simulation results discussed hereafter 

'Corresponds to the time interval in wliich tlie channel characteristics do not suffer expressive variations. 
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were obtained under static channels condition, situation where the channel gain' coefficients hold constant 
during all convergence period, i.e., for N iterations executed on each Ts\ot seconds. It is straightforward 
to show that those analyses and results can be applied considering the dynamic channels condition, where 
the channel coefficients changing following the coherence time of the channel, been observed the bound: 
Tsiot ^ {^t)c ~ -f-^ — , for all adopted mobilities, bounded by the maximal Doppler frequency fo^^^- 

TABLE I 

MULI IRAI H DS/CDMA S^ S I BM PARAMH I HRS 

Parameters Adopted Values 



DS/CDMA Power-Rate Allocation System 



Noise Power 


P„ = -63 [dBm] 


Chip rate 


Rc = 3.84 X 10^ 


Min. Signal-noise ratio 


SNR^in = 4 dB 


Max. power per user 


fma^ = 20 [dBm] 


Min. Power per user 


f min = SNR^in + P„ [dBm] 


Time slot duration 


Tsiot = 666.7ms 


# mobile terminals 


K e {5; 30} 


# base station 


BS = 4 


cell geometry 


rectangular, with Xctn = 2/ceU = 5 Km 


mobile term, distrib. 


~ U[Xcelh J/ceu] 


Channel Gain 


path loss 




shadowing 


uncorrelated log-normal, cr^ = 6 dB 


fading 


Rice: [0.6; 0.4] 


Max. Doppler freq. 


/Cmax = 11-1 Hz 


Error estimates 


G = {1+ s)G, where s ~ U[±S] 




6 = 0: 0.02 : 0.2 


User Types 


# user classes 


L = 3 (voice, data, video) 


User classes Rates 


Rmin — [^28'! 32' IG-l'^^ [I'PS] 


Verhulst Power-Rate algorithm 


Type 


partially distributed 


a 


range [0.10; 0.95] 


Optimization window 


N G [100; 1000] iterations 


Performance parameters 


Trials number, TR 


100 samples 
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Multi-cell: four base and K= 25 mobile stations location 




500 1000 1500 2000 2500 3000 3500 4000 4500 5000 
distance [m] 



Fig. 1. K = 25 Mobile terminals and 4 base stations location over 25 Km multicell rectangular geometry. 



A. Typical Convergence Performance 

Typical convergence behavior for two fixed a (slow and fast convergence scenarios), K = 7 multirate 
users, with rate assignment uniformly distributed over the three rates, [30; 120; 240] [Kbps], with = 
3.84 X 10^ chips per sec, is shown in Fig. [2l for CIR (superior plots) and power allocation solution (inferior 
plots). Plateaux indicate convergence to the optimum power vector, p*; hence dot lines (Popt in legend) 
indicates analytical solution given by dH). Note the fast convergence for all users when a = 0.9, i.e., 
^ 25 iterations, against ^ 150 for a = 0.1. Evidently, the quality of solution in this two situation is 
distinct, as discussed in the next subsection. 



Verhulst.K= /users, a = 0.1 Verhulst. K= 7users, a = 0.9 ^ .|gi Verhulst K= 7users;a = 0.1 ^ i Verhulst K= 7users; a = 0.9 




50 100 150 50 100 150 50 100 150 50 100 150 

Iteration, N Iteration, N Iteration, N Iteration, N 



Fig. 2. CIR and power convergence for K = 7; Rk.min = [120, 120, 240, 120, 30, 120, 30][Kbps]. a = 0.1 and 0.9 



B. Performance under Channel Gain Error Estimates 

In a real scenario, the SINR estimations at BS are not perfect, in the sense that the obtained values 
by estimation possess a random error characteristic. In order to incorporate this characteristic, a random 
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error is added in each element of channel gain matrix, in each iteration basis. The ratio of the estimated 
and real channel gain values is given by ^ = {l + e)gij, where e will be considered as a random variable 
with uniform distribution in the range [—5; 6]. In the subsequent simulations the adopted range values for 
6 were to 0.2, in steps of 0.02. 

1) Dependence of Solution Quality in terms of a: Since we have some idea how fast the Verhulst 
algorithm reaches the equilibrium with different values of a, it is important to determine solution quality 
in terms of convergence time. For the same system configuration of Fig. [21 we have obtained in Fig. [3] 
the associated NSE ratio, defined as: 



One can see from Fig. |3] that regardless of channel error estimates 5, the quality of solution for both 
a = 0.9 (fast) and a = 0.1 (slow convergence) at the final section of iterations (A^ > 170) shows 
high similarity {NSER ^ 1), but with a slight advantage in terms of convergence for a = 0.1. In that 
region, with both convergence factors, the algorithm approaches to the optimal solution at same speed; 
as a consequence the NSER — )• 1. Conversely, after a initial approaching convergence, i.e., after 23 and 
until ^ 120 iterations, the Verhulst algorithm with a = 0.9 produces a much better solution, resulting in 
NSE(a = 0.9) << NSE{a = 0.1). Due to the insufficient number of iterations, the algorithm is not 
able to achieve convergence for a = 0.1. 

In conclusion, the best choice for a depends on the number of iterations constraint. If the number of 
iteration is a concern, the natural choice consists in to adjust the convergence factor as high as possible. 
Otherwise low values for a produce NSE slightly smaller. 



Fig. 3. NSE Ratio considering fast (a = 0.9) and slow (a — 0.1) power convergence behavior of K = 7 mt, and 5 channel gain estimation 
error values. Bottom graph is a zoom in over A'^ G [21; 150] initial iterations. 



NSER 



NSE{a = 0.9) _ NSEifast converg.) 
NSE{a = 0.1) ~ NSE{s\ow converg.) 



(17) 



LU 
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2) Solution Quality as a Function of System Loading: Fig. |4] shows the average NSE behavior when 
the channel gain error 6 increases for the 1000th iteration and increasing system loading, K = 10, 20 and 
30 mobile terminals with different user class rates realizations (and uniformly distributed over the three 
user class rates). The convergence factor was assumed fixed a = 0.2 and the algorithm convergence ran 
100 times to each combination of K, S and user rates. 

Note from Fig. |4] that the NSE values increase for low system loading (small K), showing an increasing 
degradation rate under specific system loading when the channel error estimates S increase. We can explain 
this dependence by granularity effect, i.e., under high system loading, the average norm distances between 
the proposed algorithm solution and the analytical optimum solution results smaller due to large number 
of active users (high granularity), when compared to the low loading system cases (small K and low 
granularity). Anyway, in terms of NSE, the channel gain error estimates has a progressive effect over the 
solution quality. 



G = G(1 +e), e uniformly distrib. [-8, 5]. Avg. over TR = 1 00 samples.a= 0.2 




Fig. 4. NSE degradation as a function of K mobile terminals and S channel gain estimation error, a = 0.2. 

C. Adaptive Convergence Methods Performance 

In order to speed up the algorithm convergence, we have suggested in section Ull-Al two adaptive criteria 
based on SNIR's difference. There are two important performance aspects to be analyzed, considering 
adaptive methods against fixed a optimization methods: convergence time and solution quality. However, 
in order to privilege the quality solution analysis, in this subsection the lower convergence factor was 
adopted for the three methods; hence, in the next we evaluate the results just in terms of convergence 
time (number of iterations, A^), considering fixed the convergence factor, a = 0.1. 

A first approach to evaluate the reduction in the convergence time with the adoption of adaptive 
convergence factor is provided in Figure [51 Both plots were generated under the same channel conditions 
with = 30 users and no channel gain error estimates. As expected, the convergence rate, mainly at 
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the beginning (early iterations), is greatly speeded up. One can see that roughly the adaptive tanh— a 
procedure allows the proposed algorithm to achieve total convergence 50% early regards to the fixed 
factor a = 0.1. 



xiO"" Fixeda; K=30;a = 0;S = 



X 10" 



Tanh-a; K= 30; a = 0; 5 = 




200 400 600 800 200 400 600 800 

Iteration, N Iteration, N 

Fig. 5. Convergence speed: adaptive a with tanh function (right) and the classical fixed a method (left). 



In order to quantify the a-adaptation effect over the normalized square error, Fig. [6] shows the NSE 
and NSER for each number of iterations in the range [0;700], and considering the same scenario of 
Fig. [51 Note that for any iteration after the initial iterations (A^ > 100) the adaptive convergence factor 
provides at least one and half order better performance in terms of NSE. Specifically, for N > 200 results 
A^S'-E'(aAdpt) ~ 2- 10""^ NSE{a = 0.1). In other words, the proposed tanh a adaptive method can achieve 
the same solution quality of a; = 0.1 using less iterations 140 less iterations when > 100). 



NSE and NSER for a=0.1 and Adaptative a (Hyperbolic Tangent) K= 30 users 




100 200 300 400 500 600 700 
iteration, N 



10" 



100 200 300 400 500 600 700 
iteration, N 



Fig. 6. NSE and NSER to the adaptive method using tanh function against fixed a = 0.1. 



Finally, in order to determine the best adaptive convergence factor method suggested by (fT4l) and (fT5l) 
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we can compare the NSE for both methods under the same channel and system conditions. Fig. |7] shows 
the simulation results considering K = 30 users and 6 = 0. Note that the hyperbolic tangent mapping 
always results a lower NSE, although this difference is marginal. Therefore, for any number of iterations 
the convergence solution provided by the tanh method is better than that provided by the SNR to SNR 
target difference mapping of (fT4l) . at cost of a marginal increment in computational effort spent with tanh 
evaluation. 

SNR based and Hyperbolic Tangent Adaptativea 30 users 

10° 

10"' 
10-^ 

LU 

</) 

10-^ 

10"^ 



50 100 150 200 250 300 350 400 

Iteration, N 

Fig. 7. NSE for both a— adaptive methods. K — 30 users 

D. Computational Complexity 

As in a distributed method, each link mobile terminal-BS performs separately their updating, i.e., as 
a whole the power control is performed by K processors in parallel and each one performs only scalar 
operations. Hence, the analysis is reduced to the study of one iteration on each mobile terminal. So, 
comparing the proposed and the classical Foschini [21 algorithms, both result in same complexity. 

On the other hand, in order to compare the computational complexity of the proposed algorithm with 
analytical matrix inversion approach, we have to quantify the number of additions and multiplications 
as a function of the number of interfering mobile terminals (K — 1). Equations (fTTI) . (fT2l) and (fTSi) are 
evaluated at each iteration in each terminal using the proposed algorithm with adaptive a-tanh method. 
Table HI] shows the number of additions and multiplications operations executed per iteration. The K tanh 
evaluations per iteration were admitted as a look-up table operations. 

Hence, computational complexity of the proposed algorithm is J\f{K + 10), where M is the number 
of iterations necessary for convergence. Comparing with the best case complexity of the matrix inversion 
operation, which is given by OlK"^ ■\og{K)] [fTOll . [[TTIl . the proposed optimization methodology achieves a 
considerable complexity reduction when the number of mobile terminals is large and the NSE requirement 




13 



TABLE II 

Operations per iteration for the proposed algorithm executed on each mobile terminal. 



Equation 


Operation 


Number of Operations 


(fTTli 

^ i- i-V 


AHHiti ons 


2 




Multiplications 


3 


(Ell 


Additions 


K 




Multiplications 


K + 3 




Additions 


1 




look-up table 


1 



is not excessively tight. Besides, in the proposed method, the complexity could be controlled simply 
specifying the maximal admissible NSE. 

V. Conclusions 

An extension on discrete Verhulst power equilibrium approach, previously suggested in literature was 
proposed, taking into account the jointly power-rate optimal allocation problem. For this purpose, multirate 
users associated to voice, data and video types of traffic were aggregated as distinct user' classes, with 
the assurance of QoS and minimum rate allocation per user. Furthermore, two criteria for convergence 
speed up were suggested and compared with the fixed convergence factor case. 

Numerical results for convergence time (number of iterations), quality of solution (NSE) and number of 
basic operations (multiplications and sums) point out advantages of the Verhulst power-rate a— adaptive 
algorithm when compared to analytical solution based on the interference matrix inverse. 

Finally, the logistic map approach applied to resource allocation problem in DS/CDMA systems sug- 
gested here demonstrates tremendous potential of applicability. Future directions include a) power-rate 
allocation for multiple-input-multiple-output (MIMO) CDMA systems, and b) discrete Verhulst equilibrium 
adaptation to jointly minimize power consumption and maximize the throughput by changing the last 
constraint in (QOl) to > i?^;^. 
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